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PATENT APPLICATION 

CLASSIFYING CELLS BASED ON INFORMATION CONTAINED IN 

CELL IMAGES 

CROSS-REFERENCE TO RELATED PATENT APPLICATIONS 

5 This application claims priority from the following co-pending U.S. patent 

applications, U.S. Patent AppUcation number 09/310,879 by Vaisberg et al., and titled 
DATABASE METHOD FOR PREDICTIVE CELLULAR BIOINFORMATICS; U.S. 
Patent Application number 09/311,996 by Vaisberg et al,, and titled DATABASE SYSTEM 
INCLUDING COMPUTER FOR PREDICTIVE CELLULAR BIOINFORMATICS; and 

10 U.S. Patent Application number 09/311,890 by Vaisberg et al., and titled DATABASE 
SYSTEM FOR PREDICTIVE CELLULAR BIOINFORMATICS. Each of ttiese 
applications was filed on May 14, 1999 and is incorporated herein by reference for all 
purposes. 

BACKGROUND OF THE INVENTION 

15 The present invention pertains to image analysis methods used to classify cells based 

upon their state in the cell growth and division cycle. 

Many interesting biological conditions affect or are affected by changes in ceU cycle 
status. For example, some particular condition may cause cells to divide less JBrequently than 
normal. Another condition may cause cells to reside in the DNA synthesis state for longer 
20 than a normal amount of time. Biological "conditions" of interest to researchers include 
disease states, normal unperturbed states, quiescent states, states induced by exogenous 
agents, etc. Valuable insight may be gained by inducing a biological condition through a 
genetic manipulation, exposure to a particular agent (e.g., a compound, radiation, a field, 
etc.), deprivation of required substance, and other perturbations. 

25 In drug discovery work, valuable information can be obtained by understanding how 

a potential therapeutic affects cell growth and division. Often, this information gives some 
indication of the mechanism of action associated with the compound. For example, a 
particular class of drugs or genetic manipulations may arrest cell growth at the G2 stage 
(second gap phase). The drugs of this class are known to act via a particular set of 

30 mechanisms of action. Another class of drugs or genetic manipulations arrests cells while in 
mitosis and acts via a different mechanisms. The abihty to quickly determine whether a 
population of cells is arrested in G2 or mitosis (or some other stage) provides a valuable tool 
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in assessing the mechanism of action of an uncharacterized compound that has be&n tested 
on the population of cells. 

Conmionly, the stage of a given cell in the cell growth and division cycle is 
determined by measuring the quantity of DNA in the cell. Most cell components are made 
continuously throughout the so-called "interphase" period, between ceU divisions. However, 
DNA synthesis is an exception. DNA in the cell nucleus is repUcated oiily during a limited 
portion of the interphase, deemed the "S" phase of the cell cycle (for "synthesis")- The other 
distinct stage of the cycle is the cell-division phase, which includes both nuclear division 
(mitosis) and the cytoplasmic division (cytokinesis) that follows. The entire cell-division 
phase is denoted as the **M" phase (for ^'mitosis**)- This leaves the period between the M 
phase and the start of DNA synthesis (the S phase), which is called the "Gi" phase (first gap 
phase), and the period between the completion of DNA synthesis and the next M phase, 
which is called the "G2" phase. Interphase is thus composed of sequential Gi, S, and G2 
phases, and can comprise 90% or more of the total cell cycle time. 

Figure IB is a simple graph depicting how the quantity of nuclear DNA in a cell 
nucleus changes with cell cycle. As shown, when the cell cycle begins at the Gi phase, the 
total quantity of nuclear DNA in the cell has a value 2N. That quantity remains constant for 
the duration of Gi. At the onset of the S phase, the total quantity of nuclear DNA begins to 
increase and steadily grows. By the end of the S phase (beginning of the G2 phase), the total 
quantity of nuclear DNA has reached a value of 4N. This quantity remains constant 
throughout the G2 phase and mitosis, until two daughter cells are formed. 

Today, ceU-cycle analyses are commonly performed using fluorescence-activated 
cell analysis. This process employs a machine in which a cell suspension is forced through a 
fine nozzle and an optical measurement is made and recorded for each individual cell as it 
briefly passes through a window. Initially, a growing cell population is treated with a 
fixative (to arrest cell division and make the membranes permeable) and contacting them 
with a dye that becomes fluorescent only when it binds to DNA. When a cell is treated in 
this way, the intensity with which it fluoresces is approximately proportional to the amount 
of DNA that it contains. By passing such cells through a fluorescence analyzer, one can 
rapidly determine the relative fluorescence of the large number of cells, and, therefore, their 
relative amounts of DNA. Those cells with the least amount of DNA are in the Gi phase, 
those with double this amount are in the G2 or M phase, while cells in the S phase have 
intermediate amoxmts. The lengths of the Gi, G2 plus M, and S phases of the cell cycle can 
be calculated from the firaction of cells in each of these categories. This process is described 
in detail in various sources including Alberts et al. "Molecular Biology of the Cell" Garland 
Publishing, Inc. 1983. 
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Unfortunately, information provided by the fluorescence analyzer is too coarse for 
many applications. Most importautly, the analyzer only reads out total DNA content per 
cell; it is unable to distinguish between mitotic cells and G2 phase cells. Therefore, it would 
be desirable to have an improved process for classifying cells in a manner that distinguishes 
the interphase and mitotic states. 

SUMMARY OF THE INVENTION 

The present invention addresses this need by providing methods and apparatus that 
analyze images of cells and categorize the cells in particular ceU cycle phases based upon 
certain features extracted from the images. Independently of this, the invention can quantify 
the total amount of DNA in a cell based on specific features of an image such the as 
fluorescence intensity generated firom fluorescent molecules that bind to DNA. The 
invention can also characterize a cell as mitotic or interphase based on morphological and 
textual parameters such as the variance of the pixel intensities observed ia a cell image 
and/or the size of a region containing DNA. Of particular interest, the image analyss 
methods of this invention can accurately classify a cell into one of the following five phases: 
Gi, S, G2, early mitotic stages (prophase, prometaphase, metaphse), and late mitotic stages 
(anaphase and telophase). 

One aspect of this invention provides a method of estimating an amount of DNA in a 
cell. The method may be characterized by the following sequence: (a) receiving an image of 
the cell, in which local values of a DNA image parameter correspond to amounts of DNA at 
the locations within the cell shown on the image; and (b) estimating a total value of the DNA 
image parameter taken over a region of the cell where DNA is deemed to be present. 

The DNA image parameter indicates the amount or concentration of DNA at any 
given location in the image. In one example, the cell is treated with an agent that selectively 
associates with DNA and emits a signal recorded as the DNA image parameter. An example 
is a stain or other material that selectively binds to DNA. The DNA image parameter may 
be a Ught or radiation intensity, such as an electromagnetic radiation intensity provided at a 
particular wavelength or range of wavelengths. 

Estimating the total value of the DNA image parameter may be accompUshed by 
various techniques such as summing a per pixel value of the DNA image parameter over all 
pixels in the region of the cell where DNA is deemed to be present. In one example, the per 
pixel value of the DNA image parameter comprises a corrected per pixel intensity value such 
as the difference of an absolute pixel intensity value and an image background level. Note 
that the estimated total value of the DNA image parameter should correspond to the actual 
amount of DNA in the cell under consideration. Thus, the estimated total value of the DNA 
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image parameter may be converted to an estimated value of the total DNA in the region 
where DNA is deemed to be present. 

The method may also classify the cell into a cell cycle state based on the estimated 
value of total DNA. Classification may employ any of a number of available techniques 
5 such as applying a mixture model to operate on estimated values of total DNA for a 
population of cells. 

Another aspect of this invention pertains to an image analysis apparatus for 
estimating the amount of DNA in a cell. The apparatus may be characterized by the 
following features: (a) a memory or buffer adapted to store, at least temporarily, an image of 

10 the cell, in which image local values of a DNA image parameter correspond to amounts of 
DNA at the locations within the cell shown on the image; and (b) a processor configured or 
designed to estimate a total value of the DNA image parameter taken over at least a region of 
the cell where DNA is deemed to be present Preferably, the apparatus also includes an 
interface adapted to receive the image of the cell. In some embodiments, the system also 

15 includes an image acquisition system that produces the image of the cell. The DNA image 
parameter may be obtained and analyzed as described above. The analysis is generally 
perforaied by the processor in conjimction with the memory. 

Another aspect of the invention provides a method of distinguishing between mitotic 
and interphase cells. The method may be characterized by the following sequence: (a) 

20 receiving an image of a cell; (b) firom the image, extracting values of one or more mitosis 
indicator parameters that correspond to a cell division state of the cell; and (c) classifying the 
cell as either mitotic or interphase based upon the extracted values of the one or more 
mitosis indicator parameters. To obtain the necessary mitosis indicator parameters, the cell 
should be treated with an agent that selectively associates with DNA and emits a signal 

25 recorded as a location of DNA within the cell. Stains and other materials of the type 
described above may be employed for this purpose. 

Various types of mitosis indicator parameters will work. Examples include a 
variance in DNA concentration within the cell, the size of a region of DNA within the cell, 
and a maximal concentration of DNA within the cell. The mitosis indicator parameter or 
30 parameters is/are chosen to indicate the degree to which DNA within the cell has condensed 
into chromosomes or the degree to which DNA within the cell has separated. DNA that has 
concentrated into one or more discrete locations suggests that mitosis is occmring. A 
specific set of parameters for this purpose include (1) the average DNA concentration values 
and (2) the area occupied by DNA. 
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Various suitable techniques for classifying the cell as either mitotic or interphase 
exist. One preferred approach employs a mixture model to operate on the values of the one 
or more mitosis indicator parameters for a population of cells. In a refinement of this 
method, the algorithm classifies a mitotic cell as pre or post-anaphase. 

5 Another aspect of this iavention pertains to an image analysis apparatus for 

distinguishing between mitotic and mterphase cells. This apparatus may be characterized by 
the following features: (a) a memory or buffer adapted to store, at least temporarily, an 
image of a cell; and (b) a processor configured or designed to analyze said image in manner 
that extracts values of one or more mitosis indicator parameters that correspond to a cell 

[0 division state and classifies the cell as either mitotic or interphase based upon the extracted 
values of the one or more mitosis indicator parameters. Preferably, the apparatus includes an 
interface adapted to receive the image of the cell, and it may also include an image 
acquisition system that produces the image of the cell. The mitosis indicator parameters 
discussed above may be used by the processor to analyze and classify cell images. 

15 Yet another aspect of the invention pertains to a method of classifying a ceE based on 

an image of that cell. This method may be characterized by the following sequence: (a) 
receiving an image of a cell showing locations where nuclear DNA exists within the cell and 
also showing local amounts of the nuclear DNA; (b) fi*om the image, extracting values of 
one or more mitosis indicator parameters that correspond to a cell division state of the cell; 

20 (c) fi-om the image, estimating a total amount of DNA in the cell; and (d) classifying the ceU 
into a cell cycle phase based upon the one or more mitosis indicator parameters and total 
amount of DNA. This technique can generally classify a cell into any of the Gi, S, G2, and 

mitotic cell cycle phases. In preferred embodiments, the technique can more specifically 
classify a cell into any the Gi, S, G2, pre-anaphase mitotic, and post-anaphase mitotic 

25 phases. 

Typically, the operation of classifying the cell involves comparing the amount of 
DNA and/or the one or more mitosis indicator parameters to a model, such as a mixture 
model, providing boundaries between certain cell cycle phases in parameter space. To this 
end, the technique may also generate the model from a population of cell images. From 
30 these images, the one or more mitosis indicator parameters are extracted and the amount of 
DNA is estimated. The population used to generate the model may be an aggregation of cell 
images taken from multiple sources such as multiple wells on an assay plate. 

Still another aspect of the invention pertains to an image analysis apparatus for 
distinguishing between mitotic and interphase cells. This apparatus may be characterized by 
35 the following features: (a) a memory or buffer adapted to store, at least temporarily, an 
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image of a cell showing locations where nuclear DNA exists within the cell and also 
showing local amounts of the nuclear DNA; and (b) a processor conjBgured or designed to (i) 
analyze said image in manner that extracts values of one or more mitosis indicator 
parameters that correspond to a cell division state, (ii) estimate a total amount of DNA in the 
cell, and (iii) classifying the cell into a cell cycle phase based upon the one or more mitosis 
indicator parameters and total amount of DNA. 

It is often important to employ a sufficiently large and diverse population of cell 
images to effectively develop a classification/discrimination model. Aspects of the 
invention address this concern. Specifically, the invention provides a suitable method for 
determining parameters for classifying cells into stages of the cell cycle. This method 
involves the following operations: (a) obtaining cell images for a population of cells; (b) 
organizing the cell images based upon at least one of (i) an amount of DNA in the cells and 
(ii) a mitotic discriminator for the cell images; and (c) firom the resulting organization of cell 
images, determining the parameters for classifying the cells into stages of the ceU cycle. In 
some instances, the population of cells includes cells treated under control conditions. In 
other instances, the population of cells includes a first set of cells treated with a fiurst 
concentration of a biologically active agent and a second set of cells treated with a second 
concentration of the biologically active agent. This is refenred to as a dilution group. 
Sometimes the popxilation of cells includes both a control group and a dilution group. In a 
preiferred embodiment, detennining the parameters for classifying the cells into stages of the 
cell cycle involves applying a mixture model to the organization of cell images using an 
expectation maximization fitting technique. 

Another aspect of the invention pertains to computer program products including a 
machine readable medium on which is stored program instructions for implementing any of 
the methods described above. Any of the methods of this invention may be represented as 
program instructions that can be provided on such computer readable media. 

These and other features and advantages of the present invention will be described 
below in more detail with reference to the associated drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 A depicts the sequential phases of the cell cycle. 

Figure IB is a plot showing how the quantity of nuclear DNA varies as a function of 
cell cycle phase. 
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Figure 2 is a process flow diagram depicting - at a high level - one process of this 
invention for using an image of a cell to classify the cell according to its cell cycle phase and 
distinguishing between mitotic and interphase cells. 

Figure 3 illustrates how an image of multiple cells may be segmented to provide , 
separate representations of individual cells and thereby allow a cell-by-cell analysis. 

Figure 4 a plot of intensity versus pixel location in an image and showing how 
thresholding may be used to segment an image into individual cells. 

Figure 5 is a process flow diagram illustrating an image analysis process of this 
invention that estimates the amount of DNA in a cell using an image of the cell. 

Figure 6 is a histogram of an image used to illustrate how the 'Tjackground" of an 
image may be identified for the piupose of background subtraction. 

Figure 7 is a cartoon diagram depicting the various mitotic phases of a typical 
eucaryotic cell. 

Figure 8 is a process flow diagram illustrating an image analysis process of this 
invention that classifies cells as either mitotic or interphase cells. 

Figure 9 is a diagram showing the intensity distribution of pixels in a mitotic cell and 
in an interphase cell. 

Figure 1 0 is process flow diagram illustrating an image analysis process of this 
invention that classifies cells as one of Gi, S, G2, pre-anaphase mitotic, or post-anaphase 

mitotic. 

Figure llA is a distribution of a population of cells divided into three groups based 
on amount of nuclear DNA: (a) Gi and post-anaphase, (b) S, and <c) G2 and pre-anaphase. 

Figure 1 IB is a distribution of a population of cells mapped according to a '^mitotic 
discriminator" parameter that segregates cells into mitotic cells and interphase cells. 

Figure 12 is a distribution of a population of cells divided into five groups based 
upon amount of DNA and intensity variance: (a) Gi, (b) S, (c) G2, (d) post-anaphase, and 

(e) pre-anaphase. 

Figure 13 is a block diagram of a computer system that may be used to implement 
various aspects of this invention such as the various image analysis algorithms of this 
invention. 
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Figure 14A is a total intensity histogram showing how cells of a population are 
distributed differently when exposed to differing concentrations of Taxol. 

Figure 14B is a histogram showing how cells of a population are distributed 
differently between the mitotic and interphase stages when exposed to differing 
concentrations of Taxol. 

Figure 14C presents graphs showing how cell populations are distributed between 
Gj, S, G2, telophase, and other mitotic phases as a function of Taxol concentration. 

Figure 14D is a total intensity histogram showing how cells of a population are 
distributed when treated only with DMSO. 

Figure 14E is a histogram showing how cells of a population are distributed between 
the mitotic and interphase stages when exposed to only DMSO. 

Figure 14F presents graphs showing how cell populations are distributed between 
Gi, S, G2, telophase, and other mitotic phases within eight different wells, all treated with 

DMSO (no Taxol). 

Figure 14G is a plot showing how the cells treated with Taxol distribute themselves 
between cell cycle stages when plotted according to a mitotic discriminator and total 
intensity. 

DETAILED DESCRIPTION OF THE INVENTION 
OVERVIEW 

Generally, this invention relates to image analysis processes (methods) and apparatus 
allowing image analysis. It also relates to machine-readable media on which is provided 
instructions, data structures, etc. for performing the processes of this invention. In 
accordance with this invention, images of cells are manipulated and analyzed in certain ways 
to extract relevant features. Using those features, the apparatus and processes of this 
invention, can automatically draw certain conclusions about tlie biology of a cell. For 
example, the invention can identify the amount of DNA in a cell from its image. The 
invention can also determine which phase of the cell cycle that the cell was in when its 
image was taken. 

Relevant features from an image include morphological, textual and intensity 
features . For cell cycle analyses, relevant features often pertain to the amount and/or 
distribution of DNA in a cell, particularly nuclear DNA. The invention is not limited to 
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these features, as other cellular markers such as cytoskeletal proteins (e.g., tubulin), proteins 
having levels of expression that change significantly when the cell approaches mitosis (e.g., 
cyclin B), proteins mdergoing additional post-translational modification during mitosis (e.g,, 
histones, which are phosphorylated), certain kinesins, and certain cellular organelles (e.g., 
5 centrosomes), for example will also mark changes at certain stages of the cell division cycle. 
Thus, some embodiments of this invention consida: the amount and/or distribution of these 
other markers. Some embodiments employ a combination of features firom multiple 
markers, including DNA for example. 

The invention will now be described in terms of particular specific embodiments as 
10 depicted in the drawings. However, as will be apparent to those skilled in the art, the present 
invention may be practiced without the employing some of the specific details disclosed 
herein. Some operations or features may be dispensed with. And ofi:en alternate elements or 
processes may be substituted. 

A high level process flow 201 in accordance with one embodiment of this invention 
15 is depicted in Figure 2. As shown, the process begins at 203 where an image analysis tool 
(typically logic implemented in hardware and/or software) obtains an image showing the 
DNA firom one or more cells. Typically, an image v^U be taken from an assay plate or other 
cell support mechanism in which multiple cells are growing or stored. The image is taken in 
a manner that allows the DNA within the cells to be identified within the inxage. Preferably, 
20 the amount of DNA as fimction of position v^thin each cell is easily determined from the 
image. 

The relevant image provided at 203 is captured by an image acquisition system. In 
one embodiment, the image acquisition system is directly coupled with the image analysis 
tool of this invention. Ahematively, the image imder consideration may be provided by a 
25 remote system unaffiliated with the image acquisition system. For example, the images may 
be acquired by a remote image analysis tool and stored in a database or other rqjository until 
they are ready for use by an image analysis tool of this invention. 

After the image has been obtained at 203, the image analysis tool segments the image 
into discrete cell representations at 205, The goal of segmentation is to perform image 

30 analysis on a cell-by-cell basis. In a preferred embodiment, each representation includes 
only those pixels where the DNA of a single cell is deemed to be present. Thus, each 
representation is a bounded collection of pixels, each having an associated intensity value. 
Note that the intensity value of each pixel represents the amount of DNA at the 
corresponding location. The shape of each representation resulting from segmentation 

35 represents the boundaries within which the cellular DNA lies. In interphase cells, the DNA 
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is normally contained entirely within the nucleus. In mitotic cells, the DNA does not reside 
within a nucleus. 

Each cell (DNA) representation obtained by segmentation is separately analyzed to 
extract varioxis relevant parameters. See 207. As mentioned, most of the relevant 
parameters can be classified as either geometric parameters or statistical parameters. 
Various examples will be presented and described below. Generally these parameters are 
chosen to indicate the amoimt of DNA in a given cell and/or the distribution of that DNA at 
locations within the cell. 

After the relevant parameters have been extracted at 207, certain of these parameters 
are used to estimate a quantity of DNA and certain other parameters associated with DNA 
condensation are used to estimate a degree of condensation. See blocks 209 and 211, 
respectively. Finally, at 213, the image analysis tool uses the estimated values obtained at 
209 and 211 to classify the cells into at least one of five classifications. These classifications 
each represent a unique intersection of interphase states (Gi versus S versus G2) and division 
states (mitotic versus interphase). 

Additional operations may be performed prior to, during, or after the operations 
described Figure 2. For example, "quality control algorithms" may be employed to discard 
image data based on, for example, poor exposure, focus failures, foreign objects, and other 
imaging failures. Generally, problem images can be identified by abnormal intensities 
and/or spatial statistics. 

In a specific embodiment, a correction algorithm may be applied prior to 
segmentation to correct for changing light conditions, positions of wells, etc. In one 
example, a noise reduction technique such as median filtering is employed. Then a 
correction for spatial differences in intensity may be employed. In one example, the spatial 
correction comprises a separate model for each image (or group of images). These models 
may be generated by separately siunming or averaging all pixel values in the x-direction for 
each value of y and then separately summing or averaging all pixel values in the y direction 
for each value of x. In this manner, a paraboUc set of correction values is generated for the 
image or images under consideration. Applying the correction values to the image adjusts 
for optical system non-linearities, mis-positioning of wells during imaging, etc. 

The production of the images includes cell plating, drug dilution, drag addition and 
imaging focusing. Failures in any these systems can be detected by a variety of methods. 
For example, cell plating could fail because of a clogged tip in a delivery pipette. Such 
failxu-e can be identified by adding a fluorescent dye or bead to the cell suspension. The 
fluorescence of this dye or bead is chosen to be at a different channel (wavelength) than 

10 
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DNA the marker for the nuclei. Another potential failure could occur during drug deUvery. 
To detect such failures, one can add a fluorescent dye or bead in the compound plate before 
drug dilution. The amoimt of fluorescent dye or bead is proportional to the amount of drug. 
Yet another potential problem occurs when the focus of the image acquisition system 
5 changes during imaging. To account for such spatial biases, one can employ control wells 
containing, for example, cells with no or neutral compounds interspersed throughout the 
plate. Still another problem results from foreign objects (e.g., small dust particles) in the 
well. This can be addressed with image segmentation and statistical outher identification 
techniques. 

10 Generally the images used as the starting point for the methods of this invention are 

obtained from cells that have been specially treated and^or imaged under conditions that 
contrast the cell's DNA from other cellular components and the background of the image, hi 
the preferred embodiment, the cells are fixed and then treated with a material that binds to 
DNA and shows up in an image. Preferably, the chosen agent specifically binds to DNA, 

15 but not to most other cellular biomolecules. Preferably, the chosen agent binds 
indiscriminately with DNA, regardless of its location in the cell. In other embodiments, the 
chosen agent binds to specific sequences of DNA or regions of a chromosome. The agent 
should provide a strong contrast to other features in a given image. To this end, the agent 
should be luminescent, radioactive, fluorescent, etc. Various stains and fluorescent 

20 compounds may serve this purpose. Examples of such compounds include fluorescent DNA 
intercalators and fluorescently labeled antibodies to DNA or other nuclear component. 
Examples of fluorescent DNA intercalators include DAPI and Hoechst 33341 available from 
Molecular Probes, hic. of Eugene, Oregon. The antibodies may be fluorescently labeled 
either directly or indirectly. 

25 Various techniques for preparing and imaging appropriately treated cells are 

described in U.S. Patent Applications 09/310,879, 09/311,996, and 09/311,890, previously 
incorporated by reference. In the case of cells treated with DAPI or other fluorescent 
material, a collection of such cells is illuminated with light at an excitation frequency. A 
detector is tuned to collect Hght at an emission frequency. The collected light is used to 

30 generate the image and highlights regions of high DNA concentration. 

One approach to segmentation is depicted in Figure 3. As shown there, an image 301 
includes a pluraUty of cell images 303, each having the cell DNA highhghted. Segmentation 
effectively converts image 301 into discrete images/representations for the DNA of each cell 
as shown at 305. In a preferred embodiment, this collection of representation 305 is 
35 provided as a mask providing intensity as a fimction of position for each cell In image 301. 
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Individual cell representations 305 may be extracted from image 301 by various 
image analysis procedures. Preferred approaches include edge finding routines and 
threshold routines. Some edge finding algorithms identify pixels at locations where 
intensity is varying rapidly. For many applications of interest here, pixels contained within 
5 the edges will have a higher intensity than pixels outside the edges. Threshold algorithms 
convert all pixels below a particular intensity value to zero intensity in an image subregion 
(or the entire image, depending upon the specific algorithm). The threshold value is chosen 
to discriminate between cell (DNA) images and background. All pixels with intensity values 
above threshold in a given neighborhood are deemed to belong to a particular cell. 

10 The concepts underlying thresholding are well knovra. The technique is exemplified 

in Figure 4, which presents a plot 401 of intensity versus pixel location for an entire image 
such as image 301. For simpUcity, pixels firom a single row of an image are considered. A 
threshold value 403 is chosen to extract those features of the image having intensity values 
deemed to correspond to actual cells. In this example, peaks 405, 407, and 409 all contain 

15 collections of pixels having intensity values above threshold 403. Therefore, each of these is 
deemed to be a separate "cell" for extraction during segmentation. Because peak 411 Ues 
entirely below thresh hold 403, it is not identified as a discrete cell. 

An appropriate threshold may be calculated by various techniques. In a specific 
embodiment, the threshold value is chosen as the mode (highest value) of a contrast 

20 histogram. In this technique, a contrast is computed for every pixel in the image. The 
contrast may be the intensity difference between a pixel and its neighbors. Next, for each 
intensity value (0-255 in an eight byte image), the average contrast is computed. The 
contrast histogram provides average contrast as a function of intensity. The threshold is 
chosen as the intensity value having the largest contrast. See *The Image Processing 

25 Handbook," Third Edition, John C. Russ 1999 CRC Press LLC IEEE Press, and "A Survey 
of Thresholdmg Techniques," P.K. Sahoo, S. Soltani and A.K.C. Wong, Computer Vision, 
Graphics, and Image Processing 41, 233-260 (1988), both of which are mcorporated herein 
by reference for all purposes. 

In a specific embodiment, edge detection may involve convolving images with the 
30 Laplacian of a Guassian filter. The zero-crossings are detected as edge points. The edge 
points are Unked to form closed contours, thereby segmenting the relevant image objects. 
See The Image Processing Handbook, referenced above. 

At some point, an image analysis process must obtain image parameters relevant to 
the biological condition imder consideration. Typically, the parameters of interest relate to 
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the size, shape, contour, and/or intensity of the cell images. Examples of specific parameters 
for analyzing the cell cycle stages include the following: 



Total Litensity simi of pixel intensities in an object 

Average Intensity average intensities in an object 

Area number of pixels in an object 

Axes Ratio ratio of lengths of axes of a fitted eUipse 

Eccentricity distance from the center of an ellipse to its focus 

Solidity measure of pixels inside versus pixels outside an object 

surroimded by a simple shape 

Extent the area of the object divided by area of the smallest box to 

contain the object. 

Q05 5% quantile of total intensity 

Q25 25% quantile of total intensity 

Q50 50% quantile of total intensity 

Q75 75% quantile of total intensity 

Q95 95% quantile of total intensity 



Image analysis routines for extracting these various parameters can be designed using 
well known principles. See The Image Processing Handbook, referenced above. In 
addition, various commercially available tools provide suitable extraction routines. 
Examples of some of these products include the MetaMorph Imaging System, provided by 
Universal Imagiag Corporation, a company with headquarters in West Chester, PA and NIH 
Image, provided by Scion Corporation, a company with headquarters in Frederick, 
Maryland. 

IMAGE ANALYSIS TO ESTIMATE QUANTITIES OF DNA 

One feature of the present invention is the ability to accurately measure the amount 
of DNA in a given cell from a particular image of that cell. While the ability to accomplish 
this has multiple uses, the process depicted in Figure 2 provides a particularly useful 
apphcation. See block 209. In general, an agent appUed to highlight DNA should emit a 
signal that is proportional to the amount of agent that has bound to DNA. Thus, the amoimt 
of signal (usually indicated by the signal intensity) provides a direct indication of the 
quantity of DNA present in a given cell. When this is the case, the present invention allows 
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one to obtain an accurate measure of the amount of DNA in a given cell. Each of the above- 
mentioned stains is suitable for identifying the amount of DNA in a cell. 

Turning now to Figure 5, a specific process sequence 501 that may be used to 
measure an amount of DNA in a cell is depicted. In the context of process 201, sequence 
501 corresponds to blocks 207 and 209, As shown in Figure 5, the process begins at 503 
when the image analysis process receives a cell image. Often that image will have been 
derived fi-om a larger, segmented image of multiple cells. As indicated, the image will 
highlight the regions where DNA is present. To facilitate image analysis, the cell image is 
typically a digital representation of the cell, or a component of the cell such as its nucleus or 
other region where DNA is present. 

Next, the image analysis process determines local values of a *T)NA image 
parameter.'* See 505. Commonly the DNA image parameter will be an intensity value of 
light or radiation shown in the image. Often, the intensity value will be provided on a per 
pixel basis. In addition, the intensity value may be provided at a particular wavelength or 
narrow range of wavelengths that correspond to the emission firequency of an agent that 
specifically associates with cellular DNA. 

After the process has determined the local values of the DNA image parameter, it 
next estimates a total value of the DNA image parameter over a region of interest. See 507. 
Note that the region of interest may be an entire cell, or some subsection of the cell such as 
an organelle within the cell (e.g. the nucleus or mitochondria). While the local values of the 
DNA image parameter may correspond to local concentrations of DNA within the region of 
interest, the total value of the image parameter corresponds to the total quantity of DNA in 
the region of interest. 

In one embodiment, the process measures the total value of the DNA image 
parameter by summing a per pixel value of the DNA image parameter over all pixels in the 
region of interest. Often the per pixel values will require some correction prior to 
summation. That correction may involve, for example, subtracting a backgroxmd value 
and/or correcting for a non-linearity in the image acquisition system/method. 

After the process has estimated the total value of the DNA image parameter, it next 

converts that estimated total value to an estimated total amount of DNA in the region of 

interest. See 509. Commonly the signal presented in the image (e.g., intensity) at any given 

location is directly proportional to the amoimt of DNA at that location in the imaged target 

cell. In such cases, the conversion of block 509 involves applying a calibration fimction; in 

the simple case merely multiplying the total value of the DNA image parameter by a 

multipUer obtained through calibration. Typically, the multiplier is chosen so that the 
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amount of DNA is provided in a convenient unit such as picograms per cell. Note that 
because many applications for this aspect of the invention, require only a relative estimate of 
the amount of DNA in a region of interest, the units can be arbitrary. In cases where only a 
relative estimate of the amount of DNA is required, the process may dispense with operation 
5 509 entirely, and merely use the value obtained at operation 507, 

Finally, the analysis process outputs the value of total DNA calculated at 509. See 
511. Because this value of DNA quantity may have intrinsic value to a researcher, the 
output may simply involve displaying or printing out the value so that it can be read and 
interpreted by a researcher. Altematively, the value may be passed to a higher-level image 
10 analysis process such as process 201 depicted in Figure 2. If this is the case, the value may 
be used by another module or algorithmic operation such as that depicted in block 213 of 
process 201 . 

As mentioned, one convenient technique for determining the total amount of DNA in 
a given cell involves siunming the intensities of aU pixels provided in a given cell image. 

15 The summation is commonly referred to as the "total intensity" associated with the image. 
Other techniques may be employed to assess the amount of DNA in a given cell based on an 
analysis of that cell's image. For example, it may be useful to calculate an average intensity 
of the pixels in a cell and then scale that value for the area of the cell image. Other 
techniques for determining the amoxmt of DNA from an analysis of a cell image will be 

20 readily apparent to those you have skill in the art. 

Because any given image acquisition system will exhibit at least some variabihty and 
some non-linearity, it will often be necessary to "correct" the intensity values provided in an 
image. The correction function should generally be a monotonically non-decreasing 
function of pixel intensity. TypicaDy, it will be necessary to subtract a background value 
25 from each pixel in a ceU image. In this case, the total intensity may be given by the 
following expression: 

it pixe ls 

TOTAL INTENSITY = Y.^I,.-B) 

In this expression, the variable li represents the total intensity of any given pixel in 
the image (or a monotonically non-decreasing function of pixel intensity). The constant B 
30 represents the calculated backgroimd that must be subtracted from all pixel intensity values. 
Note that the value of B is calculated for an entire image, which often includes many 
different cells. After such image is segmented as described above, each resulting cell image 
from the parent image is analyzed using the same value of B. 
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Another way to view the total intensity of a given cellular image involves calculating 
the area under a plot of intensity versus pixel ntimber for a given cellular image. For 
example, if the level 403 shown in Figure 4 represents a background level (B) as opposed to 
a thresh hold value, then shaded areas under peeks 405, 407, and 409 represent the total 
intensity (and hence the amoimt of DNA) in each of the cells associated with those peaks. 

In one approach, a backgroimd level for a given image can be identified by analyzing 
a histogram 610 as depicted in Figure 6. As shown there, the number of pixels having a 
particular intensity is plotted versus intensity per pixel. A typical histogram will include two 
or more peaks. Usually, there will be a lower intensity and relatively large background peak. 
Smaller, higher intensity, peaks will also be discemable. In the example depicted in Figure 
6, histogram 601 includes a large, low intensity, peak 605 associated with the background of 
the image and a smaller, high intensity, peak 603 associated with the cells (DNA) of the 
image. The background level, B, subtracted from the per-pixel intensity values may be 
selected to be the intensity at the peak of curve 605. Other procediures for determining the 
value of B from a histogram such as histogram 601 will be readily apparent to tiiose who 
have skill in the art. 

Regarding other "corrections" used in measuring intensity, the absolute magnitude of 
intensity can vary from image to image due to non-linearities in the image acquisition 
procedure and/or apparatus. Specific optical aberrations can be mtroduced by various image 
collection components such as lenses, filters, beam splitters, polarizers, etc. Other non- 
linearities may be introduced by an excitation light source, a broad band Hght source for 
optical microscopy, a detector*s detection characteristics, etc. For example, some optical 
elements do not provide a "flat field." As a result, pixels near the center of the image have 
their intensities exaggerated in comparison to pixels at the edges of the image. A correction 
algorithm may be applied to compensate for this effect. Such algorithms can be easily 
developed for particular optical systems and parameter sets employed using those imaging 
systems. One simply needs to know the response of the systems under a given set of 
acquisition parameters. 

DISTINGUISHING MITOTIC AND INTERPHASE CELLS 

As indicated, it is often important to distinguish between mitotic and interphase cells. 
Depending upon how far a particular cell has progressed into mitosis, an image of the cell 
may possess a single region of DNA or two separated regions of DNA. If it possesses a 
single region, that region should have a quantity of DNA associated with the G2 phase. If it 
possesses two separate regions, each of those regions should have a quantity of DNA 
corresponding to the Gi phase. Sometimes a segmentation operation will produce two 
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separate "objects" for a late phase mitotic cell: one image for each of the two separated DNA 
regions. Thus, a simple estimation of the amount of DNA in an image segment may not be 
able to distinguish between late phase mitotic cells (in which separate regions of DNA are 
sometimes improperly identified as separate cell) and Gi interphase cells. Similarly, such 
estimation may not alone distinguish between early stage mitotic cells and interphase G2 
cells. In process 201, these distinctions are accompUshed at 211 by quantifying one or more 
parameters indicating a degree to which the DNA in the cell has condensed and/or separated. 

The significance of DNA condensation and separation can be better understood with 
reference to Figure 7, which illustrates in cartoon fashion five separate phases of mitosis. As 
a cell enters mitosis, its DNA begins to condense within the nucleus so that one can see 
individual chromosomes. This phase is referred to as the prophase of mitosis and is 
schematically illustrated by cell 707 including a nucleus 704 in which individual 
chromosomes 706 have become visible. At the end of prophase, the cell's nuclear envelope 
begins to breakdown so that the individual chromosomes now reside within the cell's 
cytoplasm. At this point, the mitotic cell has entered the prometaphase in which the 
chromosomes 706 of cell 702 have moved into the cell's cytoplasm 708. Note that the 
chromosomes are still localized near the center of cell 702. 

Next, the chromosomes within the cytoplasm tightly ahgn themselves with one 
another as shown in the third illustration of the cartoon sequence. This is the mitotic cell's 
metaphase. In the next stage in the progression (represented in the fourth image of the 
sequence), the chromosomes separate into two groups 71 OA and 71 OB while moving away 
from each other within the cell's cytoplasm 708. This is the mitotic cell's anaphase. 

Finally, the separated groups of chromosomes partition themselves into two new 
nuclei 7 12 A and 712B and the cytoplasm begins to pinch of in the middle to form two 
daughter cells. Within new nuclei 712A and 712B the chromosomes begin to de-condense 
into more evenly distributed nucleic acid. This final stage of the mitotic cell is the 
telophase. The daughter cells generated firom the telophase mitotic cell each temporarily 
reside in the Gi phase. Then they typically progress through the S and G2 phases before 
entering the next mitotic phase. 

One aspect of the present invention provides an image analysis process that can 
distinguish between interphase and mitotic cells. In a particularly preferred embodiment, the 
process distinguishes between interphase and mitotic cells by considering parameters that 
indicate the degree to which nuclear DNA has condensed to chromosomes. Examples of 
such parameters include parameters that relate to the average intensity of DNA within a cell, 
the variance in the intensity of DNA within a cell, the size of a region where the DNA 
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resides within a cell, and/or a ma?dmal concentration of DNA within a cell (often indicated 
by a maximal sigoal intensity detected in an image of the cell). 

One process, 801, for classifying cells as either mitotic or interphase is depicted in 
Figure 8. As shown there, the image analysis process begins at 804 upon receipt of a cell 
5 image showing locations of DNA, and possibly the intensity (concentration) associated with 
that DNA. Such images may be obtained via the segmentation processes described above, 
for example. 

Next, the process extracts values of one or more parameters selected to indicate a cell 
division state (e.g., either mitotic or interphase). See 806. Thus, the parameter should have 
10 a comiection to the physical location and/or concentration of DNA within a cell as a 
signature of a particular cell division state. 

As indicated in the discussion of Figure 7, mitotic cells generally have DNA that has 
condensed into chromosomes. As a consequence, the size of the region occupied by cellular 
DNA may decrease. Further, the solidity of the region occupied by the DNA may decrease 

15 due to the irregular shape of the chromosomes. Still further, the maximum intensity 
(corresponding to maximmn DNA concentration) may increase in cells due to concentration 
of DNA into discrete chromosomes. Finally, such cells may exhibit a relatively wide 
distribution in the intensity of DNA, due to the localization of DNA into the discrete 
chromosomes. Thus, examples of parameters useful to distinguish mitotic from interphase 

20 cells include size parameters (e.g., area, perimeter, diameter, length, etc.), parameters 
indicating a contour (e.g., soUdity), the average or maximum local intensity, and statistical 
parameters indicating a statistical distribution in local intensity values (e.g., the total mean, 
kurtosis, variance or standard deviation in intensity). 

The particular parameter or group of parameters used for this purpose should be 
25 chosen to work at a given set of image acquisition conditions. For example, the size and 
intensity (variance, average, or maximum) parameters work well at magnifications of about 
5x for example. At higher magnifications, sohdity can serve as a usefiil parameter for 
distinguishing between mitotic and interphase cells (by discriminating between condensed 
and uncondensed chromosomes). 

30 In one approach to discriminating between mitotic and interphase cells, the invention 

employs a mitotic discriminator that accounts for two or more relevant parameters. One 
specific approach includes both the average intensity and the area associated with the DNA. 
For example, the mitotic discriminator may be a direct function of average intensity and an 
inverse function of area (e.g., log(average intensity) minus log(area)). Cells with relatively 
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high values of this discriminator would be deemed mitotic and those with relatively low 
values of the discriminator would be deemed interphase. 

As mentioned, variance may also be a useful mitotic discriminator. The variance in 
DNA intensity may be calculated using the basic statistical expression for variance: 

5 

For a sample defined by the collection of pixels in a cell image, n is the number of 
pixels in the image, I{ is the intensity of pixel i, and I bar is the mean of the sample. Figure 9 

presents simple histograms for two hypothetical cells, one mitotic and the other interphase. 

10 As shown, the mitotic cell has a significantly wider distribution of intensities and hence a 
larger variance. The concept of variance and mitosis can be further understood by referring 
to Figure 4. As shown there, the pixels in a cell image associated with region 407 have a 
relatively wide variation in intensity (note the rough profile of its peaks) and the pixels in an 
image associated with region 405 have a relatively small variance. The individual peaks 

15 within region 407 are intended to correspond to individual chromosomes. Thus, a process of 
this invention will more likely classify the cell associated with region 407 as mitotic and the 
cell associated with region 405 as interphase. 

Returning to Figure 8, after the relevant parameter values have beeu extracted at 806, 
the image analysis process (running on an image analysis tool) next classifies the cell under 

20 consideration as either mitotic or interphase based upon the parameter values. See 808. 
Generally the cell may be classified by considering the relative magnitude of the appropriate 
values extracted at 806. The magnitude of iutensity variance is a simple example of this. 
Note however that particular parameter values may not imiversally associate with a 
particular cell division state. For example, different image acquisition techniques may 

25 introduce variability in the parameter values for similarly classified cells. Further, different 
cell types may have different parameter signatures for the same cell division state. And even 
further, cells of the same type in the same cell division state may exhibit different parameter 
values depending upon how they were grown and treated. 

Thus, it will often be desirable to generate a model from a population of similarly 
30 situated cells. Similarly situated cells generally, though not necessarily, are of the same cell 
type (e.g., the same genotype). They may also have been exposed to the same environmental 
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influences (e.g., growth media, exposure to treatment agents, physical handling, etc.). When 
such population of similarly situated cells is organized according to the values of the one or 
more parameters of interest, the image analysis process can determine how the cells 
distribute themselves. Depending upon the parameter \mder consideration, the cells should 
5 distribute themselves into at least two distinct groups: mitotic and interphase. The process 
may employ a classification technique to divide the cells of the popxdation into a Class A 
(mitotic) and a Class B (interphase). Various classification techniques are available for this 
purpose. These techniques will be discussed in more detail below. One example of such 
technique employs a mixture model. 

10 The classification operation 808 can also consider the amount of DNA in a particular 

image segment. Using this information together with the extracted parameter values 
obtained at 806, a model can fiirther classify a mitotic cell into either a pre-anaphase or a 
post-anaphase state. This process will be described in more detail below. 

Returning again to Figure 8, image analysis process 801 concludes at 810 with the 
15 process outputting the cell classification generated at 808. As with process 501, this output 
can represent the end of a complete analysis or the passage of control to another operation in 
a more involved image analysis algorithm. In the context of Figure 2, for example, the 
output provided at 810 may pass to process 201 for use in operation 213. Note further that 
operations 808 and 810 of process 801 may be performed within operation 211. Still fiirther, 
20 operations 804 and 806 of process 801 may be performed within operation 207 of higher 
level process 201. 

IMAGE ANALYSIS FOR HIGH RESOLUTION PHASE CLASSIFICATION 

As mentioned, one aspect of this invention has as its goal classifying cells into one of 
the following five states: Gl, S, G2, pre-anaphase mitotic, and post-anaphase mitotic. One 
25 implementation of this aspect of the invention is depicted in Figure 10. As shown there, an 
image analysis process 1001 performs the desired classification. 

The process begins at 1003 with the process identifying a population of similarly 
situated cells. As discussed above, these are cells that are similar in terms of their cell type 
and or their exposure to environmental influences. As such, the population of these similarly 
30 situated cells should distribute itself into various phases of the cell growth and ceU division 
cycle. In one example, the cells of interest are liver cells from a single organism that have 
been exposed to a particular drug under consideration. Further, each of these cells has been 
imaged by the same imaging system. After the relevant population of cells has been 
identified at 1003, the process next receives images of the individual cells in the population. 
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See 1005. These images may be provided as the output of a segmentation process as 
illustrated at block 205, for example. 

With the cell images now available for analysis, the image analysis process 
categorizes the cells based upon the amount of DNA that each contains. This allows the 
5 process to distinguish between those cells in the interphase states Gi, S, and G2, and 

between the mitotic states pre-anaphase and post-anaphase. See 1007. The amount of DNA 
in each cell may be estimated using the technique described in Figure 5, for example. Note 
that operation 1007 does not necessarily allow distinction between pre-anaphase and G2 or 

between post-anaphase and Gi . 

10 Preferably, the classiJBcation performed at 1007 involves a model that evaluates the 

relative positions of all cells of the population as a function of their DNA content. Assuming 
that the population of cells behaves roughly as excepted, a distribution of cells according to 
their amount of DNA should show two large groupings: one associated with GI, piost- 
anaphase, and the other associated witti G2, pre-anaphase. The cells between these two 

15 groupings will be in the S phase. 

Figure llA shows a histogram 1103 (solid line) of the number of cells in the 
population versus the amount of DNA in each cell (represented as total DNA intensity). As 
shown, the histogram 1103 includes a first peak 1105 having a relatively low amount of 
DNA and a second peak 1107 having a relatively high amount of DNA. The cells that 
20 comprise peak 1105 are either in the GI or post-anaphase state. The cells in the second 
peak, 1107, are in either the G2 or pre-anaphase state. Cells between these two peaks are in 
the S phase as indicated at saddle 1109. 

Using distributions such as that shown in Figures 11, image analysis process 1001 
(Figure 10) can classify cells based on the amount of DNA using conventional techniques. 
25 In one embodiment, operation 1007 involves the use of a mixture model, which technique is 
described in more detail below. 

After the process 1001 has performed operation 1007, it next categorizes the cells of 
the population based upon certain cell division parameters in order to distinguish between 
mitotic and interphase cells. Again, this process may involve use of a mixture model. 
30 Further, appropriate examples of process details have been set forth in the description of 
Figure 8. 

At this point in the process, the cells have been classified into growth states based 
upon the amount of DNA and have been classified into cell division states based upon a 
mitotic indicator parameter or parameters. Using this information, image analysis process 
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1001 next classifies the cells as either interphase Gi, interphase S, interphase G2, mitotic pre- 
anaphase, or mitotic post-anaphase. See 1011. 

Figure 12 shows a typical distribution of a cell population based on (a) the amount of 
DNA and (b) the size and intensity of nuclei in the cells. As shown, the ceUs can be 

5 classified into five separate groups based upon their positions in this two-dimensional 
representation. Cells with relatively low amounts of DNA and little variance in the intensity 
of their DNA label (and/or small nuclei) are classified as Gi phase cells. Cells having a 
similar degree of variance in their DNA label intensity but slightly more total DNA are 
classified as S phase cells. Cells with even more total DNA and comparable degrees of 

10 variance in the DNA label intensity are classified as G2 phase cells. Cells having relatively 
high degrees of variance in the intensity of their DNA label (and/or small nuclei) are 
classified as either telophase or pre-anaphase (including anaphase) based upon how much 
DNA they contain. Note that, as is typical, most cells cluster into the Gi phase and the G2 
phase. Note that the mitotic cells are grouped into two categories, while the interphase cells 

15 are grouped into three categories. The choice of a dividing line between pre and post- 
anaphase mitotic cells can be drawn in various ways. In one approach, the mitotic cells 
having intermediate amounts of DNA (e,g., S phase quantities of DNA) are grouped with the 
pre-anaphase mitotics. 

After the cells from the population are appropriately classified at 1011, the process 
20 outputs this classification, or a relevant component thereof, at 1013. The classification is 
output in any convenient format and may be provided to a researcher or another module in a 
higher level image analysis algorithm. In a preferred embodiment, the entire population of 
cells or some subset of those cells is used in a higher level clustering or classification 
analysis. See 1015. For example, the cell state classifications may be employed to cluster a 
25 particular gene or agent (e.g., a potential therapeutic) based on mechanism of action, 
toxicity, etc. 

Note that in the above example, the high resolution classification into one of five 
phases involved a separate classifications based upon amount of DNA and one or more other 
parameters indicative of nuclear DNA condensation. Other techniques employing other 

30 parameters for classification may be employed. For example, mitotic cells may be classified 
as early stage or late stage based upon the local distribution of DNA. If a mitotic cell's 
nuclear DNA has separated into two distinct regions, then the cell can be classified as late 
stage mitotic. If the cell's nuclear DNA has not separated into two distinct regions, then it 
can be classified as early stage mitotic. To this end, an image analysis process of this 

35 invention will extract parameters indicating the relative position of two distinct regions of 
DNA. If two regions, each having a GI quantity of DNA are observed within a defined 
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distance, this may identify a late phase mitotic cell. If a single region having a G2 quantity 
of DNA is observed, then this may suggest an early phase mitotic cell. By considering other 
parameters in addition to the relative closeness of similar regioixs of DNA, one can obtain 
even greater resolution in a classification analysis. For example, the process may be able to 
distinguish among three or more of the mitotic phases. 

CLASSIFICATION USING POPULATION STATISTICS 

In many cases, as mentioned, a given cell image can be meaningfully classified only 
when a population of similarly situated cells is considered. The members of such population 
should distribute themselves in a way that conforms with a particular distribution model. 
For example, at any given instant in time most interphase ceUs are expected to exist in the 
Gj and G2 phases, because these phases last longer than the S phase. 

When an image analysis process has properly extracted parameter values firom a 
population of cell images, it may employ any of various techniques to estimate the 
boundaries in parameter space between various classes of cells. One example of such 
technique includes the mixture model. This technique considers an expected distribution 
based upon physical model together with the actual distribution of the population. It then 
computes the estimated boundaries of parameter space. 

In a preferred embodiment of this invention, a population of ceUs is classified 
according to the quantity of DNA by applying a mixture model to various images. This may 
be appUed at 1007 in Figure 10, for example. The mixture model assumes that a given 
distribution of data points is a mixture of two or more component populations, each of which 
is represented as a normal (Gaussian) distribution. A fitting technique such as an 
expectation maximization algorithm is employed to estimate the parameters of the 
component populations. The expectation maximization algorithm is a statistical algorithm 
that maximizes the mixture likelihood using an iterative optimization routine. It is a well- 
known algorithm that is described in Dempster, A. P., Laird, N. M., and Rubin D. B., 
'Maximum Likelihood from Incomplete Data via the EM Algorithm," Journal of ttie Royal 
Statistical Society B, 1977 [39]: 1-38, for example. See also Coleman, D. A., Dong X., 
Hardin, J., Rocke D. M., and Woodruff D. L., "Some Computational Issues of Cluster 
Analysis," Joumal of Computational Statistics and Data Analysis, 1999 July; 31[1]: 1-12. 
Both of these references are incorporated herein by reference for all purposes. Note that 
maximization of the mixture likelihood is but one of many non-hierarchical methods for 
clustering or partitioning data. Other methods such as K-means may work equally well in 
some cases. 
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Inputs to the expectation maximization algorithm are the data points (e.g,, total 
intensity of labeled DNA versus number of cells), number of component populations (e.g., 
three components: Gi, S, and G2), and optionally constraints on the mean positions of each 
component (e.g., 2N for Gi, 3N for S, and 4N for G2). The outputs of the model define the 
Gaussian distribution of each component member of the population. The input data may be 
pretreated by removing outUers for example cells having much lower DNA content than 2N 
or much higher DNA content than 4N. 

The distributions of the component populations may be defined in terms of their 
means and standard deviations for example. From this information, appropriate cutoff 
positions can be defined for the overall distribution. For example a first cutoff may indicate 
the position at which it is 50% hkely that the cell is in the Gi (or post anaphase mitotic) 

phase and 50% likely that the cell is in the S phase. . 

Considering Figure llA again, the distribution 1103 (solid line) actually comprises 
three component populations (normal distributions) indicted by dashed curves: Gi/post- 
anaphase M component 1113, S component 1115, and G2/pre-anaphase M component 1117. 
Using this assumption, expectation maximization technique estimates the mean and standard 
deviation or variance for each of the three component populations. The point at which a 
given cell is equally likely to be in the S phase or the Gi/post-anaphase M phase(s) is given 

by 1 1 19. Similarly, the point at which a given cell is equally likely to be in the S phase or 
the G9/pre-anaphase M phase(s) is given by 1 121. These points are easily identified as the 

points where the Gaussians for adjacent populations intersect. 

Typically, the DNA content of the Gi phase component is centered at 2N, the DNA 
content of the S phase component is centered at 3N, and the DNA content of the G2 phase 
component is centered at 4N. In a preferred embodiment, the DNA content is represented as 
the base 2 logarithm of total intensity. This fimction has exhibits a good separation into 
three normal distributions. Using this discriminator function, the first peak and second peaks 
are centered one-half unit apart. The third peak is centered another one-half unit to the right. 

Figure 1 IB shows appUcation of a mixture model to discriminate between mitotic 
and interphase cells. This specific technique assumes that mitotic cells have small bright 
nuclei and interphase cells have larger diffiise nuclei. Consistent with this, a "mitotic 
discriminator" is defined. Such discriminator considers both the intensity and the size of the 
objects resulting from segmentation. Objects having small areas and high average intensities 
are deemed mitotic. Objects having larger areas and lower average intensities are deemed 
interphase. 
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In a specific embodiment, the mitotic discriminator is defined by the following 

expression: 

MD = log2(avg_intensity) — log2 (area) 

where "avg_intensity is the average intensity of the object under consideration and 
''area" is the total area of the object. 

If one plots MD versus number of cells, a distribution 1149 such as that shown in 
Figure 1 IB is obtained. This distribution includes a first peak 1151 and a second peak 1 153. 
Peak 1151 represents typical mitotic phase cells, while peak 1153 represents typical 
interphase cells. Employing a mixture model as described above, one may obtain a 
component 1 157 representing the mitotic population and a component 11 59 representing the 
interphase population. The Gaussian representations of these components intersect at 1161, 
the designated cutoff point between mitotic and interphase cells. 

There is a relationship between the DNA amount and the mitotic discriminator 

fimction. Note that log(total intensity) = log(average intensity) + log(area). And MD = 

log(average intensity) — log(area). The discriminators differ by the sign on the log(area) 

component. The correlation between the discriminators is [variance(log(average intensity)) 
— variance(log(area))]/square root[variance(log(average intensity))*variance(log(area))]. 

This quantity is small when the variance(log(average intensity)) is approximately equal to 
variance(log(area)). This weak correlation implies that there is very little loss in breaking 
the classification analysis into two parts. 

The abiUty of most techniques to accurately classify cell images depends on the 
number of samples in the population. Better results are obtained with larger populations. In 
some embodiments of this invention, suitably large populations may be comprised of 
samples from various sources. While it is generally preferable to use "similarly situated" 
samples as discussed above, it may sometimes be appropriate to aggregate cells fi*om slightly 
different sources. For example, it may be desirable to aggregate cells fi'om multiple wells of 
a single assay plate, where each well was treated with a different concentration of the same 
compound. In another example, the individual wells used all employ the same treatment but 
contain different strains of the same organism. 

Note that some treatment or growth regimens may have highly impredictable effects 

on the cell division cycle. In such situations, the resulting population of cell images may not 

conform to the expected physical model of the distribution of cells in various cell cycle 

phases. For example, a particular drug, with unknoAvn effects, may arrest cell growth in the 
G2 phase, thereby greatly reducing the number of mitotic cells in comparison to the typical 
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cell population- If such population was used to develop a cell cycle phase classification for 

its members, the result could be inaccurate. Therefore, a control may be employed. The 

control will use xmtreated cells (when the effects of the treatment agent are impredictable). 

Alternatively, the control may use cells treated with an agent or combination of agents 

having known effects (e.g., DMSO or Taxol (an anti-mitotic agent)) on cell growth and 

division. For some cell lines and/or relatively smaU populations, one or more of the cell 

growth stages may naturally have relatively low abundance. For example, mitotic stages 

typically are not well represented. In such cases, the control may be chosen to increase 

abundance in the naturally underrepresented stages. The resulting population will then 

present a better sample for accurate resolution by the mixture model. Note that if there were 
only G2/G2M cells in the data, the distribution of log(total intensity) would contain only one 

peak and one would be unable to determine if the peak was the Gj/GiM, S, or G2/G2M 
peak. 

hi a particularly preferred embodiment, a calibrated cell population includes a 
dilution series. The dilution series includes multiple dilution groups, each including, for 
example, 8 wells on a plate, hi this approach, all samples include the same cell line and 
same dmg, but the drug is diluted over a wide range of concentrations. At the extreme of 
complete dilution, the cells are left untreated, except for exposure to DMSO. The cell 
images firom each level of dilution are aggregated to obtain the entire population for analysis 
using a mixture model. This ensiures a sufficiently large population for generating the 
mixture model. It increases the likelihood that all phases of the cell cycle are represented. It 
also reduces the likelihood that one or more phases firom the cell cycle will be inadequately 
represented because a particular concentration of the drug strongly suppresses one or more 
cell growth stages. Thus, the chances of a calibration/annotation error are reduced. 

A specific embodiment employs two separate control groups, one to provide 

discrimination between interphase stages and the other to provide discrimination between 

mitotic and interphase cells. The first group contains DMSO (a solvent used in all treatment 

wells and considered a neutral treatment). The second group contains Taxol (an anti-mitotic 

agent). A group dilution is run on each of tiiese control groups. In the DMSO dilution 
group, the cells will contain all three groups: Gj/GiM, S, and G2/G2M. The cutoffs 

between these groups (cutoffs A and B) can be saved from an analysis of this data. The 
other dilution group, containing Taxol, will arrest cells in mitosis at high concentrations. 
The distribution of MD over the Taxol is a mixture of two normal distributions. Running the 
expectation maximization algorithm on this data, one can obtain a cutoff between mitotic 
and interphase cells (cutoff C). Together the cutoffs from both dilution groups is used to 
classify cells. Cells with log(total intensity) < A and MD < C are classified as Gl cells. 
Cells with a log(total mtensity) < A and MD > C as classified as late phase mitotic cells. 
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Cells with A < log(total intensity) < B are classij&ed as S phase cells. Cells with log(total 
intensity) > B and MD < C are classified as G2 cells. And finally, cells with log(total 
intensity) > B and MD > C are classified as early phase mitotic cells. After all cells are so 
classified, proportions of each of the five classes are computed on a well by well basis. 

5 SOFTWARE/HARDWARE 

Generally, embodiments of the present invention employ various processes involving 
data stored in or transferred through one or more computer systems. Embodiments of the 
present invention also relate to an apparatus for performing these operations. This apparatus 
may be specially constructed for the required purposes, or it may be a general-purpose 

10 computer selectively activated or reconfigured by a computer program and/or data structure 
stored in the computer. The processes presented herein are not inherently related to any 
particular computer or other apparatus. In particular, various general-purpose machines may 
be used with programs written in accordance with the teachings herein, or it may be more 
convenient to construct a more specialized apparatus to perform the required method steps. 

15 A particular structure for a variety of these machines will appear from the description given 
below. 

In addition, embodiments of the present invention relate to computer readable media 
or computer program products that include program instructions and/or data (including data 
structures) for performing various computer-implemented operations. Examples of 

20 computer-readable media include, but are not limited to, magnetic media such as hard disks, 
floppy disks, and magnetic tape; optical media such as CD-ROM disks; magneto-optical 
media; semiconductor memory devices, and hardware devices that are specially configured 
to store and perform program instructions, such as read-only memory devices (ROM) and 
random access memory (RAM). The data and program instructions of this invention may 

25 also be embodied on a carrier wave or other transport mediima. Examples of program 
iustructions include both machine code, such as produced by a compiler, and files containing 
higher level code that may be executed by the computer using an interpreter. 

Figure 13 illustrates a typical computer system that, when appropriately configured 
or designed, can serve as an image analysis apparatus of this invention. The computer 

30 system 1300 includes any number of processors 1302 (also referred to as central processing 
units, or CPUs) that are coupled to storage devices including primary storage 1306 (typically 
a random access memory, or RAM), primary storage 1304 (typically a read only memory, or 
ROM). CPU 1302 may be of various types including microcontrollers and microprocessors 
such as programmable devices (e.g., CPLDs and FPGAs) and improgrammable devices such 

35 as gate array ASICs or general purpose microprocessors. As is well known in the art, 
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primary storage 1304 acts to transfer data and instructions uni-directionally to the CPU and 
primary storage 1306 is used typically to transfer data and instructions in a bi-directional 
manner. Both of these primary storage devices may include any suitable computer-readable 
media such as those described above. A mass storage device 1308 is also coupled bi- 

5 directionally to CPU 1302 and provides additional data storage capacity and may include 
any of the computer-readable media described above. Mass storage device 1308 may be 
used to store programs, data and the like and is typically a secondary storage medium such 
as a hard disk. It will be appreciated that the information retained within the mass storage 
device 1308, may, in appropriate cases, be incorporated in standard fashion as part of 

10 primary storage 1306 as virtual memory. A specific mass storage device such as a CD-ROM 
1314 may also pass data imi-directionally to the CPU. 

CPU 1302 is also coupled to an interface 1310 that connects to one or more 
input/output devices such as such as video monitors, track balls, mice, kej^oards, 
microphones, touch-sensitive displays, transducer card readers, magnetic or paper tape 

15 readers, tablets, styluses, voice or handwriting recognizers, or other well-known input 
devices such as, of course, other computers. Finally, CPU 1302 optionally may be coupled 
to an external device such as a database or a computer or telecommunications network using 
an extemal connection as shown generally at 1312. With such a connection, it is 
contemplated that the CPU might receive infonnation from the network, or might output 

20 information to the network in the course of performing the method steps described herein. 

hi one embodiment, the computer system 1300 is directly coupled to an image 
acquisition system such as an optical imaging system that captures images of cells. Digital 
images from the image generating system are provided via interface 1312 for image analysis 
by system 1300. Alternatively, the images processed by system 1300 are provided from an 

25 image storage source such as a database or other repository of cell images. Again, the 
images are provided via interface 1312. Once in the image analysis apparatus 1300, a 
memory device such as primary storage 1306 or mass storage 1308 buffers or stores, at least 
temporarily, digital images of the cell. Typically, the cell images will show locations where 
DNA exists within the cells. In these images, local values of a DNA image parameter (e.g., 

30 radiation intensity) correspond to amounts of DNA at the locations within the cell shown on 
the image. With this data, the image analysis apparatus 1300 can perform various image 
analysis operations such as distinguishing between mitotic and interphase cells, estimating 
the amount of DNA in a cell, and classifying a cell in a particular cell division or cell growth 
state. To this end, the processor may perform various operations on the stored digital image. 

35 For example, it may analyze said image in manner that extracts values of one or more 

mitosis indicator parameters that correspond to a cell division state and classifies the cell as 

either mitotic or interphase based upon the extracted values of the one or more mitosis 
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indicator parameters. Alternatively, or in addition, it may estimate a total value of the DNA 
image parameter taken over at least a region of the cell where DNA is deemed to be present. 

EXAMPLES 

Figures 14A-14G depict the results of an experiment in which human lung cancer 
5 epitheUal cells were treated with Taxol and the effect of Taxol on the cell cycle was 
characterized. The specific cells used in this study were fi-om Cell Line A549 (human lung 
cancer epithelial cells) (ATCC:CCL-1S5). Two day staged cell cultures of the A549 cells 
were trypsinized for five minutes from T175 cm flasks. The cells were then suspended in 
ten milUliters of RPMI media with ten percent serum and counted using both a 
10 hemocytometer and Coulter Coimter. The suspension was diluted further in media to ensure 
that there were 1600 cells per twenty microliters in suspension. The cells were then counted 
again. Thereafter, they were transferred to a Cell Stir and kept in suspension while being 
plated into barcoded 384 well plates using a MultiDrop at twenty microhters per well. The 
plates were transferred to a hmnidified carbon dioxide incubator to recover for twenty-four 
15 hours before the addition of Taxol or DMSO. Prior to imaging the cells were treated with 
DAPI to stain DNA. Cells firom each well were then imaged. 

Of relevance to this example, eight wells were treated with Taxol at varying 

concentrations and eight other "control" wells were treated with DMSO, which is not toxic 

to the cells under consideration. The concentration of Taxol varied fi-om approximately 0.09 

20 micromolar to 200 micromolar. Figure 14A shows how the concentration of Taxol affects 

the distribution of cells based upon the amount of DNA contained in each cell. Cells in each 

of the eight Taxol wells were separately imaged and the total measured intensity in each cell 

was recorded. As discussed above, the total intensity corresponds to the amount of DNA in 
each cell. A histogram of cell count versus log2 (total intensity) was prepared as shown in 

25 Figure 14A. In the histogram, the cell count has been normalized so that the area imder each 
curve was equal to a value of one. 

As shown, the Taxol concentration has a strong influence on the amount of DNA 

contained in a typical cell. At low concentrations, the cells in a sample assume a near typical 
distribution in which there are large numbers of Gi and G2 phase cells. As explained above, 

30 this condition is represented by two peaks in the total intensity histogram. The peak with the 
lower intensity represents Gj phase cells and the peak with the higher intensity represents 

G2 phase cells. Initially, at low concentrations of Taxol, the size of the G j peak is larger. 

As the concentration of Taxol increases, the second peak becomes more pronounced with 
respect to the first peak. At very high Taxol concentrations, there are very few Gj cells. In 
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other words, the relative numbers of G2 and Gj phase cells become skewed in favor G2 
phase cells. This suggests that Taxol arrests cells in the G2 phase. 

Figure 14B shows how Taxol affects a mitotic discriminator. The data used to 

produce the histogram of Figure 14B was taken from the same images of the eight Taxol 

5 containing wells used to produce the histogram of Figure 14A. In this histogram, the mitotic 
discriminator is defined as log2(average intensity) minus log2(area), A cell count histogram 

based on this mitotic discriminator typically shows two peaks, one for mitotic cells (higher 
value of MD) and another for interphase cells (lower value of MD). As shown in Figure 
14B, the higher concentrations of Taxol tends to greatly increase the relative incidence of 
10 mitotic cells. 

Figure 14C shows two graphs depicting the relative proportions of cells in each of 
five cell cycle phases as a function of Taxol concentration. The five phases are Gj, S, G2, m 

(telophase), and M (prophase, metaphase, and anaphase). The cell classifications were 
derived from the data used to construct the histograms of Figures 14A and 14B, The graph 

15 on the left used a "dilution group" to classify the individual cells of the images. The graph 
on the right used a "control" to classify the individual cells. The control employed cell 
images from the DMSO treated wells. As discussed above, a dilution group employs cell 
images from each of a series of dilutions to yield a population that is then used to find cutoff 
points between each of the various cell cycle phases. A control algorithm employs a set of 

20 control images for this purpose. The control images are derived using no agent (Taxol) 
and/or an agent having a known affect (DMSO). Note that the cell classifications developed 
using each of the two methods are remarkably similar. 

In the graphs, there is a curve with points labeled "o." These are defined as the 
number of objects divided by the average number of DMSO objects. The number of DMSO 

25 objects is the number of "cells" identified in an image of a well. This number is averaged 
over the eight DMSO wells to obtain the denominator. The numerator is the number of 
objects identified in each image of the Taxol containing wells. As shown in Figure 14C, the 
nimiber of such objects decreases with increasing concentration of Taxol. Because these 
objects correspond to living cells, the "o" curves indicate that more lung cancer cells die 

30 with increasing concentrations of Taxol, 

Figure 14D is a histogram showing cell count versus log2(total intensity) for each of 

eight images taken from eight different wells, each treated with DMSO and no Taxol. These 

images were used for the control algorithm discussed above. Note the two characteristic 
peaks for G\ and G2 phase cells. Note also the strong similarities in the curves from each of 

35 the eight images. 
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Figure 14E is a histogram showing cell count versus the mitotic discriminator 
(described above) for each of the eight images of Figure 14D. As expected, most of the cells 
reside in the interphase state and only a few reside in the mitotic state. Further, proportions 
of cells in these states are quite similar for each of the eight images. 

Figure 14F presents two graphs depicting the relative proportions of cells m each of 
five cell cycle phases as a function of well number. These graphs are similar to those shown 
in Figure 14C, but there is no change in the concentration of applied agent; all use a single 
concentration of DMSO. Not surprisingly, the relative proportions of cells in each of the 
five phases remain relatively constant. The total number of objects, o, varies slightly about 
the mean. 

Figure 14G is a scatter plot showing how image objects distribute themselves as a 
fimction of log2 (total intensity) and MD (log2(average intensity) minxis log2(area)). The 

cutoff points between the various cell cycle stages, derived as discussed above, are shown m 
the Figure. 

Although the above has generally described the present invention according to 
specific processes and apparatus, the present invention has a much broader range of 
applicabiUty. In particular, the present invention is not limited to a particular kind of data 
about a particular cell, but can be apphed to virtually any cellular data where an 
understanding about the workings of the cell is desired. Thus, in some embodiments, the 
techniques of the present invention could provide information about many different types or 
groups of cells, substances, and genetic processes of all kinds. Of course, one of ordinary 
skill in the art would recognize other variations, modifications, and alternatives. 
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CLAIMS 

what is claimed is: 

I . A method of estimating an anioimt of DNA in a cell, the method comprising: 

5 receiving an image of the cell, in which local values of a DNA image parameter 

correspond to amounts of DNA at the locations within the cell shown on the image; and 

estimating a total value of the DNA image parameter taken over at least a region of 
the cell where DNA is deemed to be present. 

10 2. The method of claim 1 , wherein the image of the cell is a digital representation of the 
cell. 

3. The method of claim 1, wherein the cell is treated with an agent that selectively 
associates with DNA and emits a signal recorded as the DNA image parameter. 

15 

4. The method of claim 3, wherein the agent is a DNA stain. 

5. The method of claim 1, wherein the DNA image parameter is a light or radiation 
intensity. 

20 

6. The method of claim 1, wherein the DNA image parameter is an electromagnetic 
radiation intensity provided at a particular wavelength or range of wavelengths. 

7. The method of claim 1, wherein estimating the total value of the DNA image 

25 parameter comprises summing a per pixel value of the DNA image parameter over all pixels 
in the region of the cell where DNA is deemed to be present. 

8. The method of claim 7, wherein the per pixel value of the DNA image parameter 
comprises a corrected per pixel intensity value. 

30 

9. The method of claim 8, wherein the corrected per pixel intensity value comprises the 
difference of an absolute pixel intensity value and an image backgroxmd level. 

10. The method of claim 9, wherein the corrected per pixel intensity value comprises the 
35 difference corrected for non-linearities in an image acquisition system used to produce the 

image of the cell. 

I I. The method of claim 1, further comprising converting the estimated total value of the 
DNA image parameter to an estimated value of the total DNA in the region where DNA is 

40 deemed to be present. 

12. The method of claim 1 1 , further comprising classifying the cell into a cell cycle state 
based on the estimated value of total DNA. 

45 13. The method of claim 12, wherein classifying the cell into a cell cycle state comprises 
using a mixture model to operate on estimated values of total DNA for a population of cells. 
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14. The method of claim 1, wherein the region of the cell where DNA is deemed to be 
present is the cell nucleus. 

15. A computer program product comprising a machine readable medium on which is 
provided program instructions for estimating an amount of DNA in a cell, the instructions 
comprising: 

code for receiving an image of the cell, in which local values of a DNA image 
parameter correspond to amounts of DNA at the locations within the cell shown on the 
image; and 

code for estimating a total value of the DNA image parameter taken over at least a 
region of the cell where DNA is deemed to be present. 

16. The computer program product of claim 15, wherein the image of the cell is a digital 
representation of the cell. 

17. The computer program product of claim 15, wherein the cell is treated with an agent 
that selectively associates with DNA and emits a signal recorded as the DNA image 
parameter. 

18. The computer program product of claim 17, wherein the agent is a DNA stain. 

19. The computer program product of claim 15, wherein the DNA image parameter is a 
Ught or radiation intensity. 

20. The computer program product of claim 15, wherein the DNA image parameter is an 
electromagnetic radiation intensity provided at a particular wavelength or range of 
wavelengths. 

21. The computer program product of claim 15, wherein estimating the total value of the 
DNA image parameter comprises summing a per pixel value of the DNA image parameter 
over all pixels in the region of the cell where DNA is deemed to be present. 

22. The computer program product of claim 21, wherein the per pixel value of the DNA 
image parameter comprises a corrected per pixel intensity value. 

23. The computer program product of claim 22, wherein the corrected per pixel intensity 
value comprises the difference of an absolute pixel intensity value and an image background 
level. 

24. The computer program product of claim 23, wherein the corrected per pixel intensity 
value comprises the difference corrected for non-linearities in an image acquisition system 
used to produce the image of the cell. 

25. The computer program product of claim 15, further comprising program instmctions 
for converting the estimated total value of the DNA image parameter to an estimated value 
of the total DNA in the region where DNA is deemed to be present. 

26. The computer program product of claim 25, further comprising classifying the cell 
into a cell cycle state based on the estimated value of total DNA. 
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27. The computer program product of claim 26, wherein classifying flie cell into a cell 
cycle state comprises using a mixture model to operate on estimated values of total DNA for 
a population of ceUs. 

28. The computer program product of claim 1, wherein the region of the cell where DNA 
is deemed to be present is the cell nucleus. 

29. An image analysis apparatus for estimating the amount of DNA in a cell, the 
apparatus comprising: 

a memory or buffer adapted to store, at least temporarily, an image of the cell, in 
which image local values of a DNA image parameter correspond to amounts of DNA at the 
locations within the cell shown on the image; and 

a processor configured or designed to estimate a total value of the DNA image 
parameter taken over at least a region of the cell where DNA is deemed to be present. 

30. The apparatus of claim 29, further comprising an interface adapted to receive the 
image of the cell. 

31 . The apparatus of claim 29, further comprising an image acquisition system that 
produces the image of the cell. 

32. The apparatus of claim 29, wherein the cell is treated with an agent that selectively 
associates with DNA and emits a signal recorded as the DNA image parameter. 

33. The apparatus of claim 29, wherein the DNA image parameter is a Ught or radiation 
intensity. 

34. The apparatus of claim 29, wherein the processor estimates the total value of the 
DNA image parameter by siunming a per pixel value of the DNA image parameter over all 
pixels in the region of the cell where DNA is deemed to be present. 

35. The apparatus of claim 34, wherein the per pixel value of the DNA image parameter 
comprises a corrected per pixel intensity value. 

36. The apparatus of claim 35, wherein the corrected per pixel intensity value comprises 
the difference of an absolute pixel intensity value and an image backgroimd level, which 
difference is corrected for non-linearities in an image acquisition system used to produce the 
image. 

37. The apparatus of claim 29, wherein the processor also converts the estimated total 
value of the DNA image parameter to an estimated value of the total DNA in the cell. 

38. A method of distinguishing between mitotic and interphase cells, the method 
comprising: 

receiving an image of a cell; 

firom the image, extracting values of one or more mitosis indicator parameters that 
correspond to a cell division state of the ceU; and 

classifying the cell as either mitotic or interphase based upon the extracted values of 
the one or more mitosis indicator parameters. 
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39. The method of claim 38, wherein the image of the cell is a digital representation ot 
the cell. 

40. The method of claim 38, wherein the cell is treated with an agent that selectively 
associates with DNA and emits a signal recorded as a location of DNA within the cell. 

41 . The method of claim 38, wherein the signal is an electromagnetic radiation intensity. 

42. The method of claim 38, wherein the agent is a DNA stain. 

43. The method of claim 38, wherein the one or more mitosis indicator parameters 
include at least one of a variance in DNA concentration within the cell, the size of a region 
occupied by DNA v/ithin the cell, an average concentration of DNA within the cell, and a 
maximal concentration of DNA within the cell. 

44. The method of claim 38, wherein the one or more mitosis indicator parameters 
include at least a variance in DNA concentration within the cell. 

45. The method of claim 38, wherein classifying the cell as either mitotic or interphase 
comprises evaluating the one or more mitosis indicator parameters to determine a degree to 
which DNA within the cell has separated. 

46. The method of claim 38, wherein classifying the cell as either mitotic or interphase 
comprises evaluating the one or more mitosis indicator parameters to determine a degree to 
which DNA within the cell has condensed into chromosomes. 

47. The method of claim 38, wherein classifying the cell as either mitotic or interphase 
comprises evaluating the one or more mitosis indicator parameters to determine a degree to 
which DNA within the cell has concentrated into one or more discrete locations. 

48. The method of claim 38, further comprising classifying a mitotic cell as pre or post- 
anaphase. 

49. The method of claim 38, wherein classifying the cell as either mitotic or interphase 
comprises using a mixture model to operate on the values of the one or more mitosis 
indicator parameters for a population of cells. 

50. The method of claim 38, wherein the image of the cell shows locations where the 
DNA exists within the cell. 

51. A computer program product comprising a machine readable medimn on which is 
provided instructions for distinguishing between mitotic and interphase cells, the instructions 
comprising: 

code for receiving an image of a cell; 

code for extracting, from the image, values of one or more mitosis indicator 
parameters that correspond to a cell division state of the cell; and 

code for classifying the cell as either mitotic or interphase based upon the extracted 
values of the one or more mitosis indicator parameters. 
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52. The computer program product of claim 5 1 , wherein the image of the cell is a digital 
representation of the cell. 

53. The computer program product of claim 51, whereiu the cell is treated with an agent 
that selectively associates with DNA and emits a signal recorded as a location of DNA 
within the cell. 

54. The computer program product of claim 5 1 , wherein the signal is an electromagnetic 
radiation intensity. 

55. The computer program product of claim 5 1 , wherein the agent is a DNA stain. 

56. The computer program product of claim 5 1 , wherein the one or more mitosis 
indicator parameters include at least one of a variance ia DNA concentration within the cell, 
the size of a region occupied by DNA within the cell, an average concentration of DNA 
within the cell, and a maximal concentration of DNA within the cell. 

57. The computer program product of claim 51 , wherein the one or more mitosis 
indicator parameters include at least a variance in DNA concentration within the cell. 

58. The computer program product of claim 5 1 , wherein classifying the cell as either 
mitotic or interphase comprises evaluating the one or more mitosis indicator parameters to 
determine a degree to which DNA within the cell has separated. 

59. The computer program product of claim 51, wherein classifying the cell as either 
mitotic or interphase comprises evaluating the one or more mitosis indicator parameters to 
determine a degree to which DNA within the cell has condensed into chromosomes. 

60. The computer program product of claim 51, wherein classifying the cell as either 
mitotic or interphase comprises evaluating the one or more mitosis indicator parameters to 
determine a degree to which DNA within the cell has concentrated into one or more discrete 
locations. 

61 . The computer program product of claim 5 1 , further comprising program instructions 
for classifying a mitotic cell as pre or post-anaphase. 

62. The computer program product of claim 51, wherein classifying the cell as either 
mitotic or interphase comprises using a mixture model to operate on the values of the one or 
more mitosis indicator parameters for a population of cells. 

63. The computer program product of claim 5 1 , wherein the image of the cell shows 
locations where the DNA exists within the cell. 

64. An image analysis apparatus for distinguishing between mitotic and interphase cells, 
the apparatus comprising: 

a memory or buffer adapted to store, at least temporarily, an image of a cell; and 
a processor configured or designed to analyze said image in manner that extracts 
values of one or more mitosis indicator parameters that correspond to a cell division state 
and classifies the cell as either mitotic or interphase based upon the extracted values of the 
one or more mitosis indicator parameters. 
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65. The apparatus of claim 64, further comprising an interface adapted to receive the 
image of the cell. 

66. The apparatus of claim 64, further comprising an image acquisition system that 
produces the image of the cell. 

67. The apparatus of claim 64, wherein the cell is treated with an agent that selectively 
associates with DNA and emits a signal recorded as a location of DNA within the cell. 

68. The apparatus of claim 64, wherein the one or more mitosis indicator parameters 
include at least one of a variance in DNA concentration within the cell, the size of a region 
occupied by DNA within the cell, an average concentration of DNA within the cell, and a 
maximal concentration of DNA within the cell. 

69. The apparatus of claim 64, wherein the one or more mitosis indicator parameters 
include at least a variance in DNA concentration within the cell. 

70. The apparatus of claim 64, wherein the processor classifies the cell as either mitotic 
or interphase by evaluating the one or more mitosis indicator parameters to determine a 
degree to which DNA within the cell has condensed into chromosomes. 

71. The apparatus of claim 64, wherein the processor classifies the cell as either mitotic 
or interphase by evaluating the one or more mitosis indicator parameters to determine a 
degree to which DNA within the cell has concentrated into one or more discrete locations. 

72. The apparatus of claim 64, wherein the image of the cell shows locations where the 
DNA exists vdthin the cell. 

73. A method of classifying a cell based on an image of that cell, the method comprising: 
receiving an image of a cell showing locations where nuclear DNA exists within the 

cell and also showing local amounts of the nuclear DNA; 

fi-om the image, extracting values of one or more mitosis indicator parameters that 

correspond to a ceU division state of the cell; 

fi-om the image, estimating a total amount of DNA in the cell; and 
classifying the cell into a cell cycle phase based upon the one or more mitosis 

indicator parameters and total amount of DNA. 

74. The method of claim 73, wherein the cell cycle phases into which the cell can be 
classified include Gj, S, G2, and mitotic. 

75. The method of claim 74, wherein the cell cycle phases into which the cell can be 
classified include Gi, S, G2, pre-anaphase mitotic, and post-anaphase mitotic. 

76. The method of claim 73, wherein classifyLug the cell comprises comparing at least 
one of the amount of DNA and the one or more mitosis indicator parameters to a model 
providing boundaries between certain cell cycle phases in parameter space. 
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77. The method of claim 76, furfher comprising generating the model from a population 
of cell images, for which the one or more mitosis indicator parameters have been extracted 
and the amoimt of DNA has been estimated. 

78. The method of claim 77, wherein the model is generated using a mixture model. 

79. The method of claim 77, wherein the population used to generate the model is an 
aggregation of cell images taken from multiple sources. 

80. The method of claim 79, wherein the multiple sources are multiple wells on an assay 
plate. 

81 . The method of claim 77, wherein the population used to generate the model is a 
collection of images of control cells having a known distribution of cell cycle phases. 

82. A computer program product of classifying a cell based on an image of that cell, the 
method comprising: 

receiving an image of a cell showing locations where nuclear DNA exists within the 
cell and also showing local amounts of the nuclear DNA; 

from the image, extracting values of one or more mitosis indicator parameters that 
correspond to a cell division state of the cell; 

from the image, estimating a total amount of DNA in the cell; and 

classifying the cell into a cell cycle phase based upon the one or more mitosis 
indicator parameters and total amount of DNA. 

83. The computer program product of claim 82, wherein the cell cycle phases into which 
the cell can be classified include Gj, S, G2, and mitotic. 

84. The computer program product of claim 83, wherein the cell cycle phases into which 
the cell can be classified include Gi, S, G2, pre-anaphase mitotic, and post-anaphase mitotic. 

85. The computer program product of claim 82, wherein classifying the cell comprises 
comparing at least one of the amount of DNA and the one or more mitosis indicator 
parameters to a model providing boundaries between certain cell cycle phases in parameter 
space. 

86. The computer program product of claim 85, further comprising program instractions 
for generating the model from a population of cell images, for which the one or more mitosis 
indicator parameters have been extracted and the amount of DNA has been estimated. 

87. The computer program product of claim 86, wherein the model is generated using a 
mixture miodel. 

88. The computer program product of claim 86, wherein the population used to generate 
the model is an aggregation of cell images taken from multiple sources. 

89. The computer program product of claim 88, wherein the multiple sources are 
multiple wells on an assay plate. 
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90. The computer program product of claim 86, wherein the population used to generate 
the model is a collection of images of control cells having a known distribution of cell cycle 
phases. 

91 . An image analysis apparatus for distmguishing between mitotic and interphase cells, 
the apparatus comprising: 

a memory or buffe adapted to store, at least temporarily, an image of a cell showing 
locations where nuclear DNA exists within the cell and also showing local amounts of the 
nuclear DNA; and 

a processor configured or designed to (a) analyze said image in manner that extracts 
values of one or more mitosis indicator parameters that correspond to a cell division state, 
(b) estimate a total amount of DNA in the cell, and (c) classifying the cell into a cell cycle 
phase based upon the one or more mitosis indicator parameters and total amount of DNA. 

92. The apparatus of claim 91, further comprising an interface adapted to receive the 
image of the cell. 

93. The apparatus of claim 91, further comprising an image acquisition system that 
produces the image of the cell. 

94. The apparatus of claim 91, wherein the cell is treated with an agent that selectively 
associates with DNA and emits a signal recorded as a location of DNA within the cell. 

95 . The apparatus of claim 9 1 , wherein the one or more mitosis indicator parameters 
include at least one of a variance in DNA concentration within the cell, the size of a region 
occupied by DNA within the cell, and a maximal concentration of DNA within the cell. 

96. The apparatus of claim 91, wherein the processor classifies the cell as either mitotic 
or interphase by evaluating the one or more mitosis indicator parameters to determine a 
degree to which DNA within the cell has condensed into chromosomes. 

97. The apparatus of claim 91, wherein the processor estimates the total amount of the 
DNA in the cell by summing a per pixel value of a DNA image parameter over all pixels in 
the region of the cell where DNA is deemed to be present, wherem the DNA image 
parameter corresponds to amounts of DNA at the locations within the cell shown on the 
image. 

98. The apparatus of claim 97, wherein the cell is treated with an agent that selectively 
associates with DNA and emits a signal recorded as the DNA image parameter. 

99. The apparatus of claim 97, wherein the per pixel value of the DNA image parameter 
comprises a corrected per pixel intensity value. 

1 00. The apparatus of claim 99, wherein the corrected per pixel mtensity value comprises 
the difference of an absolute pixel intensity value and an image background level. 

101. A computer assisted method of detemiining parameters for classifying cells into 
stages of the cell cycle, the method comprising: 

obtaining cell images for a population of cells; 
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organizing the cell images based upon at least one of (i) an amount of DNA in the 
cells and (ii) a mitotic discriminator for the cell images; and 

from the resulting organization of cell images, detemiining said parameters for 
classifying the cells into stages of the cell cycle. 

1 02. The method of claim 101, wherein the population of cells includes cells treated under 
control conditions, 

1 03 . The method of claim 101, wherein the population of cells includes a first set of cells 
treated with a first concentration of a biologically active agent and a second set of cells 
treated with a second concentration of the biologically active agent. 

1 04. The method of claim 1 03, wherein the population of cells fiorther includes a third set 
of cells treated under control conditions. 

105. The method of claim 101, wherein determining said parameters for classifying the 
cells into stages of the cell cycle comprises applying a mixture model to the organization of 
cell images. 

106. The method of claim 105, wherein determining said parameters for classifying the 
cells into stages of the cell cycle comprises employing an expectation maxindzation fitting 
technique to the mixture model. 

1 07. The method of claim 101, wherein said parameters for classifying the cells into 
stages of the cell cycle allow classification of cells into at least the following stages: Gj, S, 

G2, and mitotic. 

108. The method of claim 101, wherein said parameters for classifying the cells into 
stages of the cell cycle allow classification of cells into at least the following stages: Gj, S, 

G2, and pre-anaphase mitotic and post-anaphase mitotic. 

109. The method of claim 101, wherein the cell images are organized based upon both the 
amount of DNA in the cells and the noitotic discriminator for the cell images. 

110. The method of claim 101, wherein the mitotic discriminator includes at least one of a 
variance in DNA concentration within the cell, the size of a region occupied by DNA within 
the cell, and a maximal concentration of DNA within the cell. 

111. A computer program product comprising a machine readable medium on which is 
provided instructions for determining parameters for classifying cells into stages of the cell 
cycle, the instructions comprising: 

code for obtaining cell images for a population of cells; 

code for organizing the cell images based upon at least one of (i) an amount of DNA 
in the cells and (ii) a mitotic discriminator for the cell images; and 

code for determining said parameters for classifying the cells into stages of the cell 
cycle based on the organization of cell images. 

112. The computer program product of claim 111, wherein the population of cells includes 
cells treated xmder control conditions. 
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113. The computer program product of claim 111, wherein the population of cells includes 
a first set of cells treated with a first concentration of a biologically active agent and a 
second set of cells treated with a second concentration of the biologically active agent. 

1 14. The computer program product of claim 1 13, wherein the population of cells further 
includes a third set of cells treated under control conditions. 

115. The computer program product of claim 111, wherein the code for determining said 
parameters for classifying the cells into stages of the cell cycle comprises code for applying 
a mixture model to the organization of cell images. 

116. The computer program product of claim 115, wherein the code for detomining said 
parameters for classifying the cells into stages of the cell cycle comprises code for 
employing an expectation maximization fitting technique to the mixture model. 

117. The computer program product of claim 111, wherein said parameters for classifying 
the cells into stages of the cell cycle allow classification of cells into at least the following 
stages: Gj, S, G2, and mitotic. 

118. The computer program product of claim 111, wherein said parameters for classifying 
the cells into stages of the cell cycle allow classification of cells into at least the following 
stages: Gj , S, G2, and pre-anaphase mitotic and post-anaphase mitotic. 

119. The computer program product of claim 111, wherein the cell inoiages are organized 
based upon both the amoimt of DNA in the cells and the mitotic discriminator for the cell 
images. 

120. The computer program product of claim 111, wherein the mitotic discriminator 
includes at least one of a variance in DNA concentration within the cell, the size of a region 
occupied by DNA within the cell, and a maximal concentration of DNA within the cell. 
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